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ABSTRACT 

On the basis of the near infrared observations of bulge red clump stars near 
the Galactic center, we have determined the galactocentric distance to be R = 
7.52 ±0.10 (stat) ±0.35 (sys) kpc. We observed the red clump stars at | I |< 1?0 
and 0?7 <| b |< 1?0 with the IRSF 1.4 m telescope and the SIRIUS camera in the 
H and K$ bands. After extinction and population corrections, we obtained (m — 
M) = 14.38 ± 0.03 (stat) ±0.10 (sys). The statistical error is dominated by the 
uncertainty of the intrinsic local red clump stars' luminosity. The systematic error 
is estimated to be ±0.10 including uncertainties in extinction and population 
correction, zero-point of photometry, and the fitting of the luminosity function 
of the red clump stars. Our result, Rq = 7.52 kpc, is in excellent agreement with 
the distance determined geometrically with the star orbiting the massive black 
hole in the Galactic center. The recent result based on the spatial distribution of 
globular clusters is also consistent with our result. In addition, our study exhibits 
that the distance determination to the Galactic center with the red clump stars, 
even if the error of the population correction is taken into account, can achieve an 
uncertainty of about 5 percent, which is almost the same level as that in recent 
geometrical determinations. 
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Subject headings: Galaxy: center — Galaxy: fundamental parameters — stars: 
distances 



1. INTRODUCTION 

Since the review "The distance to the center of the Galaxy" by Reid (1993), the accuracy 
of the distance to the Galactic center (GC) from us, Rq, has been improved especially 
in the primary measurements, i.e., without a "standard candle". Salim & Gould (1999) 
showed that one can measure Rq with an accuracy of a few percent by observing the proper 
motions and radial velocities of several stars near the super massive black hole (SMBH) 
at the GC for 15—30 yr. Eisenhauer et al. (2003, 2005) actually determined Rq with a 
4 percent uncertainty by using one of the stars orbiting the SMBH. Therefore, Rq would 
be determined more precisely than that estimated from the luminosity distance, i.e., the 
secondary measurements. 

However, does the position of the SMBH coincide with the dynamical center and/or the 
luminous center? Several lines of evidence for a displacement of the centroid of the Galactic 
bar from the center of mass of our Galaxy exist (Blitz 1994; Morris & Serabyn 1996, and 
references therein), and the off-centering of the stellar bar is predicted from the simulation of 
the Milky Way model (e.g., Fux 1999) although the compact radio source Sgr A* is nearly at 
rest at the dynamical center of the Galaxy (Reid & Brunthaler 2004). Off-centered bars are 
commonly observed in external galaxies (e.g., Levine & Sparke 1998, and references therein). 
The displacement would also reveal itself as a global tilt of the inner disk of the Galaxy. A 
tilt of the triaxial bulge has been suggested (Blitz & Spergel 1991), but the studies using 
COBE data sets could not find statistically significant evidence for the tilt (e.g., Dwek et al. 
1995). The understanding of these effects is thus still in its earliest stages. 

Red clump (RC) giants have been recently claimed to be a very accurate distance 
indicator (e.g., Demers, Irwin, & Gambu 1994), and a good tracer of the Galactic bar (e.g., 
Stanek et al. 1997; Nishiyama et al. 2005). RC stars are the equivalent of the horizontal- 
branch stars for a metal-rich population, have narrow distributions in luminosity and color, 
and consequently occupy a distinct region in the color magnitude diagram (CMD). The 
Hipparcos catalog (ESA 1997) allows us an exact calibration of an RC absolute magnitude; 
therefore, RC stars can be used as a reliable standard candle. The stellar density and the 
gravitational potential of the Galaxy are well traced by low- and intermediate-mass, evolved 
stars because they are dynamically relaxed, and constitute the largest fraction of the total 
stellar mass. Thus, by comparing the distance to the central stellar cluster with that obtained 
by the distribution of bar components, the existence of the displacement can be examined. 
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From photometric observations in infrared wavelength, i? was already derived by Alves 
(2000) and Babusiaux & Gilmore (2005). However, Alves (2000) used RC stars in Baade's 
window, which is ~ 4° rs 560 pc (if R Q = 8 kpc) away from the GC, and his error in R 
is dominated by the small number (about 200) of RC stars. Babusiaux & Gilmore (2005) 
observed (I, b) = (0?0, +1?0), which is only 140 pc from the GC. Since their purpose was to 
study the bar structure, not to determine the distance to the GC, they did not provide an 
exact value of R . Therefore our determination, which employs RC stars at ] I |< 1° and 
| b |< 1° and is almost free from the statistical error due to abundant RC stars in the bulge, 
is the most reliable ever studied with RC stars. 



2. OBSERVATION AND DATA REDUCTION 

Our observations were conducted in 2002 March— July and 2003 April— August using 
the near-infrared camera SIRIUS (Simultaneous Infrared Imager for Unbiased Survey; Na- 
gashima et al. 1999; Nagayama et al. 2003) on board the IRSF (Infrared Survey Facility) 
telescope. IRSF is a 1.4 m telescope constructed and operated by Nagoya University and 
SAAO (South African Astronomical Observatory) at Sutherland, South Africa. The SIRIUS 
camera can provide J (1.25 jttm), H (1.63 /xm), and Kg (2.14 jim) images simultaneously, 
with a field of view of 7'7 x 7'.7 and a pixel scale of C'45. 

About lOOx 3 (J,H,K S ) images were obtained over | I |< 1?0 and | b |< 1?0 (Fig. 1). 
We observed only on photometric nights, and the typical seeing was l'/l FWHM in the H 
band. A single image comprises 10 dithered 5 s exposures. 

Data reduction was carried out with the IRAF (Imaging Reduction and Analysis Facil- 
ity) 1 software package. Images were prereduced following the standard procedures of near- 
infrared arrays (dark frame subtraction, flat-fielding, and sky subtraction). Photometry, 
including point-spread function (PSF) fitting, was carried out with the DAOPHOT package 
(Stetson 1987). We used the DAOFIND task to identify point sources, and the sources were 
then utilized for PSF-fitting photometry in the ALLSTAR task. About 20 sources were used 
to construct the PSF for each image. 

Each image was calibrated with the standard star 9172 (Persson et al. 1998), which 
was observed every hour in 2002 and every half- hour in 2003. We assumed that star 9172 



1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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has H = 12.12, and K s = 12.03 in the IRSF/SIRIUS system. The average of the zero- 
point uncertainties was about 0.03 mag in the two bands. The averages of the 10 a limiting 
magnitudes were H = 16.6 and K s = 15.6. 

To analyze the distribution of RC stars, we define the extinction-free magnitude 

K H _ K = K S - e{ h-K s ) X [{H ~ Ks) ~ {H ~ KsU (1) 

where we use the reddening law Ak s /E(H — Ks) = 1-44 (Nishiyama et al. 2006), and the 
intrinsic H — K$ color of RC stars (H — Ks)o = 0.07 (Bonatto, Bica, & Girardi 2004). 
Thus, K H _ K is defined so that if A Ks /E(H — K s ) is independent of location, then K H _ K is 
independent of extinction for any particular star. 

We construct Ks vs. H — Ks CMDs and extract the stars in the region of the CMDs 
dominated by RC stars. The extracted stars are used in turn to make Kh-k histograms 
(luminosity functions, see Fig. 2). The histograms have clear peaks and are fitted with 
the sum of exponential and Gaussian functions {thick curve in Fig. 2). Owing to highly 
nonuniform interstellar extinction over the region observed, the peaks of RC stars in the 
CMDs shift from one line of sight to another over the range 13.0 < K s < 14.5 although a 
dispersion of K H _ K is small (see §3.3). Since the mean H and K s magnitudes of RC stars 
become too faint in highly reddened fields, estimates of the peak magnitudes and the colors 
of RC stars can be unreliable in such fields. To avoid this problem, we use only those fields 
in which the peak magnitude of RC stars is more than 1 mag brighter than the 10a limiting 
magnitudes (Fig. 1). In addition, we confirmed the completeness to be 85% at K s = 15 by 
adding artificial stars into the most crowded image. 



3. RESULTS 

3.1. The Distance to the Galactic Center 

By fitting the luminosity function of the extracted stars at | I |< 1?0 and 0?7 <| b |< 1?0, 
we obtained the center of the RC peak as K H -k = 12.855 ± 0.005 (Fig. 2). The distance 
modulus to the GC is given by 

(m - M) = K H _ K - M Ks + AM K , (2) 

where M Ks is the absolute K s magnitude of local RC stars, and AM K is the population 
correction calculated from theoretical stellar evolution models. Here we adopt AM K = —0.07 
(scaled solar metallicity, Salaris & Girardi 2002). M Ks is estimated as M Ks = M K + 
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Fig. I. — Top left: Fields used for data analysis, plotted in Galactic coordinates. Top right: 
Distribution of the fields in Galactic latitude. Bottom left: Distribution of the fields in 
Galactic longitude. 



AMk s ~k, where Mk is the absolute K magnitude and AMk s ~k is the difference between 
magnitudes in the K s and K bands. Alves et al. (2002) obtained M K = —1.60 ± 0.03 mag, 
and Bonatto, Bica, & Girardi (2004) showed that AM Ks _ K « 0.01, and thus M Ks = -1.59. 
Hence, we obtain (m — M) = 14.38 ± 0.03 (stat), which corresponds to R = 7.52 ± 0.10 
(stat) kpc, and the resulting statistical error is the sum of the squared errors. Systematic 
errors are estimated below. 
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Fig. 2. — Histogram of the dereddened Kh-k magnitude for all the stars in the RC- 
dominated region at | I |< 1?0 and 0?7 <| b |< 1?0. Exponential and Gaussian functions are 
used to fit the histograms (thick curve). 

3.2. Error Estimation 

It is well known that the statistical error in the determination of distances using RC stars 
is very small (e.g., Paczyhski & Stanek 1998). From the fitting shown in Fig. 2, we obtain 
the center of the RC peak of Kh-k = 12.855 ± 0.005. As shown below, systematic errors 
are much larger than this, and thus our statistical error is negligible in this determination. 
However, the determination of the absolute magnitude of the local RC stars has a statistical 
error of 0.03 (Alves et al. 2002), and the total statistical error o" sta t in {jn — M)q is therefore 
estimated to be <T sta t = 0.03 mag. 

One of the systematic errors arises from the fitting of the RC peak, and thus we estimate 
the error by changing the fitting range of the exponential and Gaussian functions. The lower 
limit used in the fitting is changed from 10.5 to 12.0 in Kh-k, and the upper limit is changed 
from 13.1 to 13.5. When the upper limit is larger than 13.6, the fitting fails because of the 
slope change around 13.6 (see, Fig. 2). The smallest and largest values of the peak center 
are 12.824 and 12.885, respectively, and thus we estimate that the uncertainty in the fitting 
is ±0.03 mag. 
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In the analysis of the interstellar extinction law, Nishiyama et al. (2006) showed that 
the law could change even in the near infrared wavelengths; at | I \< 2° and | b |< 1°, 
Ak s I 1 Eh-k s changes by ~ ±0.09. Since the mean extinction correction is ~ 0.6 in H — Ks, 
the systematic error in the extinction correction is estimated to be ±0.05 mag. 

A systematic error in our observation comes from the uncertainty of the zero-point 
determination of photometry. From the analysis of the standard stars, we obtained the 
mean zero-point uncertainty to be about 0.04 mag in K H _ K . 

The population correction AM K was provided by the analysis of the RC brightness 
with the theoretical stellar evolution model (Salaris & Girardi 2002). We adopted AM K = 
^locai _ M^s e — —0.07 with solar metallicity. The systematic error for the population 
correction is difficult to estimate because the bulge star formation rate and age-metallicity 
relation should be involved. [For the Large Magellanic Cloud, Salaris, Percival, & Girardi 
(2003) estimated the systematic error.] In this study, we adopt the systematic error for the 
correction ±0.07 (Percival & Salaris 2003), which is the 1 a error of the mean residual in 
AM K between the model predicted and observed RC magnitudes. Since Salaris, Percival, 
& Girardi (2003) obtained the total systematic error in the distance modulus to the Large 
Magellanic Cloud as t^l, the systematic error ±0.07 in only AM K might be overestimated. 
Detailed analysis of the uncertainty and empirical tests will provide us with more accurate 
AM K . 

We can estimate the systematic error in system transformation by comparing isochrones 
in different filter systems. The isochrones presented in Girardi et al. (2002) are available 
online 2 . The difference in RC magnitudes and colors between the Johnson-Cousins-Glass 
system (Bessell & Brett 1988) and the 2MASS system 3 are ~ 0.01 in K$ (also mentioned 
in Bonatto, Bica, & Girardi, 2004) and ~ 0.007 in H — K$, respectively. The filter systems 
of 2MASS and SIRIUS are more similar than those of 2MASS and the Johnson- Cousins- 
Glass system, and thus the systematic error coming from the difference in the filter system 
is smaller than 0.01, which is negligible in our analysis. 

Since the fields we observed are very crowded, stars might be blended and appear 
brighter than they really are. To examine this effect, we checked for the dependence of 
the RC peaks on star density using the 38 fields shown in Fig. 1 (~ 300 and ~ 200 stars 
arcmin -2 [> 10cr, Kg band] for the most crowded and uncrowded regions, respectively), and 
found no dependence. Moreover, we confirmed that the completeness is 85% at K s = 15. 



2 http://pleiadi. pd.astro.it 

3 http://www. ipac.caltech.edu/2mass/releases/second/doc/sec3_lbl. html 
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The crowding effect is thus very small in this study. 
Table 1: Systematic Error Budget 



Error 


Estimation (mag) 


fitting the RC peak 


0.03 


extinction law 


0.05 


zero-point uncertainty 


0.04 


population correction 


0.07 


system transformation 


< 0.01 



The systematic errors in our analysis are summarized in Table 1. Combining these 
errors, we obtained the total systematic error a sys = 0.10. Therefore, we obtain (m — M)q = 
14.38 ± 0.03 (stat) ±0.10 (sys), which corresponds to #0 = 7.52 ± 0.10 (stat) ±0.35 (sys) 
kpc. 

3.3. R Dispersion 

The Galactic latitude and longitude distributions of 38 fields we used in our analysis 
are shown in Fig. 1. The averages of longitude in positive and negative Galactic latitudes 
are 0.95 (0.11 rms uncertainty) and —0.92 (0.05 rms), and those of latitude in positive and 
negative Galactic longitudes are —0.16 (0.52 rms) and 0.09 (0.58 rms), respectively. Hence, 
the fields we observed are distributed nearly symmetrically around the GC, which provides 
a good sample to examine the distribution of RC stars. 

Using the 38 fields shown in Fig. 1, we obtained the peaks of RC stars in Kh-k for each 
field, and made plots as functions of Galactic latitude {left panels in Fig. 3) and longitude 
{right panels). The top panels in Fig. 3 are the plots for the fields at b > 0°, and the bottom 
panels are at b < 0°. The straight lines in the plots are the least-squares fits to the data 
points. The slopes given by the fits are -0.04 ± 0.13 for b > 0° and -0.17 ± 0.20 for b < 0° 
as a function of Galactic latitude, and 0.02 ± 0.01 for b > 0° and 0.05 ± 0.02 for b < 0° as a 
function of Galactic longitude. These slopes suggest that the distribution of RC stars along 
the lines of sight is almost flat in our observed fields. 

The RC peak distribution for all the fields is shown in Fig. 4 {white histogram). From 
the histogram, we determine the mean K H _ K to be 12.85 and the rms uncertainty of the 
distribution to be 0.06. Because the statistical errors of the RC peak in K h ~k are very small 
(< 0.01), this uncertainty comes from the systematic uncertainties in the fitting of the K H _ K 
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distributions (0.03; see Table 1), the extinction law (0.05), and the zero-point of standard 
stars (0.04). The sum of these squared errors is 0.07, consistent with the rms uncertainty of 
0.06 obtained above. This consistency shows the validity of our estimate of the systematic 
errors. 

The dotted and hatched histograms in Fig. 4 show the RC peak distributions in K h ~k 
at b > 0° and b < 0°, respectively. The mean Kh-k an d uncertainty for b > 0° are 12.87 
and 0.05, and those for b < 0° are 12.83 and 0.07, respectively. These also indicate that the 
distances to bulge RC stars in positive and negative Galactic latitudes are almost the same. 

4. DISCUSSION 

We calibrated the absolute magnitude of RC stars by using the theoretical prediction 
for the stars in Baade's window (Salaris & Girardi 2002). Different populations could make a 
difference in the absolute magnitudes between our observational fields and Baade's window. 
However, the bulk of stellar populations in the bulge is old (e.g., van Loon et al. 2003), and 
the difference in age coming from the star formation histories employed in the models affects 
Mk much less because the RC magnitude fades very slowly for old populations (see also 
Girardi & Salaris 2001). In addition, Ramirez et al. (2000) examined a metallicity gradient 
along the minor axis of the inner bulge (0° < b < —4° at I ~ 0°) with the result that the 
gradient is negligible. We therefore conclude that the RC absolute magnitudes in our fields 
and Baade's window are the same. 

Our result R = 7.52 ± 0.10 (stat) ±0.32 (sys) kpc, is in excellent agreement with 
7.62 ± 0.32 kpc (Eisenhauer et al. 2005), which is determined geometrically with the star 
S2 orbiting the SMBH in the GC. Recently, Bica et al. (2006) provided a new estimate of 
Ro based on symmetries of the spatial distribution of 116 globular clusters, and obtained an 
average value of Rq = 7.2 ± 0.3 kpc. This is smaller than those derived by Eisenhauer et al. 
(2005) or adopted in recent reviews (8.0 ± 0.5 kpc by Reid [1993], 7.9 ± 0.2 kpc by Nikiforov 
[2004]), but is consistent with our result. In Bica et al. (2006), numerous globular clusters are 
employed, and the variation of total to selective absorption Ry was taken into account; thus 
the uncertainty in R is a factor of ~ 3 smaller than previous ones using globular clusters. 

Previously, Rq was determined from the K and K$ band photometry of RC stars (Alves 
2000; Babusiaux & Gilmore 2005; Nishiyama et al. 2005). Assuming M^ s = —1.61, Babu- 
siaux & Gilmore (2005) obtained Rq = 7.6 ± 0.15 kpc without the population correction 
(AM K ). Because AM K for the Galactic bulge makes the value of R smaller than that 
without AM K , the result of Babusiaux & Gilmore (2005) suggests that R is significantly 
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Fig. 3. — Plots of RC peaks in Kh~k as functions of Galactic latitude (left panels) and 
longitude (right panels). Top: Fields at b > 0°. Bottom: Fields at b < 0°. Straight lines are 
least-squares fits to the data points. 
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Fig. 4. — Histograms of RC peaks in Kh-k for all the fields {white histogram), the fields at 
b > 0° {dotted histogram) and those 6 < 0° {hatched histogram). The mean Kh-k and rms 
are 12.85 and 0.06 (all), 12.87 and 0.05 {b > 0°), and 12.83 and 0.07 {b < 0°), respectively. 

smaller than 8 kpc. By contrast, Alves (2000) derived R = 8.24 ± 0.42 kpc, which is ~ 2a 
larger than our result. However, he did not adopt the population effect; with AM K = —0.07, 
his distance becomes Ro = 7.98 kpc. The error in his distance modulus is dominated by 
a statistical one (±0.11 mag) because of the small number of RC stars in Baade's window. 
When the statistical error is taken into account, his distance is consistent with our result. 
Assuming {H — Ks)q = 0.1, Nishiyama et al. (2005) also derived Rq ~ 8.3 kpc without the 
population correction. When {H — Ks)o = 0.07 and AMk = —0.07 are adopted, their result 
is in good agreement with the current estimate. 

Is the center of the RC distribution more distant than 8 kpc (from the recent review by 
Reid 1993) or 8.5 kpc (1AU 1985 standard; Kerr & Lynden-Bell 1986)? When AM K = -0.07 
is applied, the distance derived by Alves (2000) and Babusiaux & Gilmore (2005) with RC 
stars also become smaller than 8 kpc. The star formation rate and age-metallicity relation 
of the Galactic bulge is extremely complicated, and theoretical prediction is troublesome. In 
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addition, metallicity dependence of AM K becomes considerable at an age larger than ~ 7 Gyr 
(Salaris & Girardi 2002), and there are few old and metal rich clusters to be compared with 
theoretical prediction (Grocholski & Sarajedini 2002; Percival & Salaris 2003). However, to 
be i?o ~ 8.0 kpc, we require AMx = +0.5, which is much larger and would be unacceptable 
from the theoretical model. Therefore, we conclude that the center of the RC distribution 
is at most nearer than 8 kpc. 

Some observations suggest a displacement of the centroid of a bar from the central 
stellar cluster in our Galaxy. A mode of (to, /) = (1,1), where a star completes / orbits while 
executing to epicycle oscillations as seen in the rotating frame, could make the displacement 
(Sellwood & Wilkinson 1993). The mode also reveals itself as a global tilt of the inner disk 
of the Galaxy. Several pieces of evidence exist for the mode of (to, Z) = (1, 1), such as the 
pronounced longitudinal asymmetry of molecular line emissions, and the tilt of the plane 
defined by gas xi orbits with respect to the large-scale Galactic plane (Blitz 1994; Morris 
& Serabyn 1996, and references therein). The studies of bulge structure (e.g., Stanek et al. 
1997; Nishiyama et al. 2005) showed that the RC star is a good tracer of the bar structure in 
the central region of our Galaxy. Thus, the mode of (to, I) = (1, 1) could make a difference 
in the distances determined from the distribution of RC stars and from the stars orbiting 
the SMBH at the center of our Galaxy. The good agreement between our result and that of 
Eisenhauer et al. (2005) indicates that the displacement along the line of sight is negligible 
with an uncertainty of a few hundred pc. 



5. CONCLUSION 

Using bulge red clump stars, the distance to the Galactic center is determined. Our 
distance modulus is estimated to be (to — M) = 14.38 ± 0.03 (stat) ±0.10 (sys), which 
corresponds to the distance R = 7.52 ± 0.10 (stat) ±0.35 (sys) kpc. The statistical error is 
dominated by that of the local red clump magnitude used for calibration. The main sources of 
our systematic error are from uncertainties in extinction correction and population correction 
calculated from theoretical stellar evolution models. Our analysis shows that the distances 
to the red clump stars are almost the same in the range of | I |< 1?0 and 0?7 <| b |< 1?0. 
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